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Quantum sensors based on solid-state spins provide tremendous opportunities in a wide range of 
fields from basic physics and chemistry to biomedical imaging. However, integrating them into a 
scanning probe microscope to enable practical, nanoscale quantum imaging is a highly challenging 
task. Recently, the use of single spins in diamond in conjunction with atomic force microscopy 
techniques has allowed significant progress towards this goal, but generalisation of this approach 
has so far been impeded by long acquisition times or by the absence of simultaneous topographic 
information. Here we report on a scanning quantum probe microscope which solves both issues, by 
employing a nano-spin ensemble hosted in a nanodiamond. This approach provides up to an order 
of magnitude gain in acquisition time, whilst preserving sub-100 nm spatial resolution both for the 
quantum sensor and topographic images. We demonstrate two applications of this microscope. We 
first image nanoscale clusters of maghemite particles through both spin resonance spectroscopy and 
spin relaxometry, under ambient conditions. Our images reveal fast magnetic field fluctuations in 
addition to a static component, indicating the presence of both superparamagnetic and ferromagnetic 
particles. We next demonstrate a new imaging modality where the nano-spin ensemble is used 
as a thermometer. We use this technique to map the photo-induced heating generated by laser 
irradiation of a single gold nanoparticle in a fluid environment. This work paves the way towards 
new applications of quantum probe microscopy such as thermal/magnetic imaging of operating 
microelectronic devices and magnetic detection of ion channels in cell membranes. 


In recent years, quantum sensing has attracted increas¬ 
ing attention because it opens up the possibilities for 
non-invasive sensing with nanoscale mapping potential 
and single atom sensitivity. The model system for these 
studies has been the negatively charged nitrogen-vacancy 
(NV) centre in diamond, which has shown great promise 
as an atomic-sized sensor operating under ambient condi¬ 
tions w- The NV centre is a photostable point defect 
whose electronic ground state is a spin triplet that can 
be initialised and read out by optical means [31 [5]. This 
makes it possible to detect the spin resonance of a single 
NV centre and coherently control its spin state. Owing 
to interaction with the local environment, measurement 
of the spectral and temporal properties of an NV cen¬ 
tre’s spin provides access to a number of local properties, 
such as the magnetic field, electric field, pressure or tem¬ 
perature [SHU]. Furthermore, NV sensors are capable 
not only of measuring static fields but also fluctuating 
ones. This, in turn, allows retrieval of spectral informa¬ 
tion from a sample P3HIS]. This capability was recently 
used to perform nuclear magnetic resonance spectroscopy 
on a small ensemble of protons pan]. 

In order to obtain spatial information about the sam¬ 
ple’s properties, the NV probe must be scanned relative 
to the sample. This is commonly achieved by using an 
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atomic force microscope (AFM), which enables a probe 
to be scanned relative to a sample, while precisely con¬ 
trolling the probe-sample distance. Integrating an NV 
probe into an AFM can be achieved by grafting a nan¬ 
odiamond onto the tip of the AFM [T9H22] . or by fab¬ 
ricating an all-diamond AFM tip hosting an NV centre 
[23l 24 . An alternative route is to make use of an NV 
centre located near the surface of a bulk diamond and to 
place the sample on the AFM tip [25H27] . Of these var¬ 
ious strategies, attachment of the nanodiamond to the 
AFM tip is the only one which can provide the detailed 
topography of the sample in addition to quantum sensor 
information. This is a key requirement for a practical, 
versatile instrument capable of imaging complex samples 
such as assemblies of nanoparticles, nanostructured thin 
films or biological cells. Furthermore, most probe mi¬ 
croscopy experiments to date have employed a single NV 
centre as the scanning probe, rather than an ensemble. 
This is motivated by the fact that the ultimate spatial 
resolution is in principle given by the size of the probe 
volume, which is minimal for a single NV centre (~ 1 
nm), although in practice the resolution is limited by the 
probe-sample distance, which is typically several tens of 
nm at least Hang. With an ensemble of NV centres, the 
spatial resolution is degraded by spatial averaging over 
the ensemble. However, the acquisition time is also dic¬ 
tated by the number N of NV centres employed in the 
sensor, since this limits the signal-to-noise ratio which 
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FIG. 1. Schematic of the scanning nano-spin ensemble microscope, (a) The probe consists of a nanodiamond containing a small 
ensemble of electronic spins (symbolised by the red arrows) grafted onto the tip of an atomic force microscope (AFM). Optical 
excitation and readout, combined with microwave (MW) excitation, enable quantum measurements of the spin ensemble 
properties such as spin resonance, relaxation or decoherence. Scanning the sample relative to the probe produces images 
of the sample with magnetic (B field), electric (E field) or temperature (T) contrast, depending on the probing technique 
used, (b) Optically detected spin resonance spectrum of a nanodiamond on a tip. The solid line is the fit to a sum of two 
Lorentzian functions centred at frequencies v± — D ± E, providing the zero-field splitting parameters D — 2868.3 ±0.1 MHz 
and E — 5.4 ±0.1 MHz. (c) Spin relaxation (left) and spin decoherence (right) curve of a nanodiamond on the AFM tip. The 
inset depicts the sequence of laser (green) and MW (blue) pulses employed. Solid lines are fits to a single exponential decay, 
revealing a spin relaxation time Ti = 142 ± 9 fis and spin coherence time T 2 = 780 ± 30 ns. 


scales as y/N [2]. As a consequence, acquisition times up 
to several tens of hours per image are typically required 
with a single NV centre, which imposes serious techni¬ 
cal constraints owing to mechanical and thermal drifts 
present in the microscope [25ti28] . 

In this work, we seek to optimise the trade-off between 
acquisition time and resolution by employing a small NV 
ensemble as our probe. Based on this idea, we demon¬ 
strate a scanning nano-spin ensemble microscope that is 
capable of providing simultaneous topographic and quan¬ 
tum sensor information with a spatial resolution better 
than 100 nm and typical acquisition times of an hour. 
We illustrate the versatility of our microscope by apply¬ 
ing various sensing modes to two different samples. We 
first use spin resonance spectroscopy and spin relaxom- 
etry to image the magnetic properties of aggregates of 
maghemite nanoparticles (Fe 2 C >3 in the cubic 7 phase). 
Secondly, we demonstrate nanoscale thermal imaging of a 
photo-heated gold nanoparticle, through spin resonance 
spectroscopy of the quantum probe in a fluid environ¬ 
ment. These results demonstrate that our approach, 
which relies on nanodiamonds containing a small ensem¬ 
ble of NV spins, is a promising avenue to realise a prac¬ 
tical, versatile scanning quantum probe microscope. 

The experimental set-up is depicted in Fig. [T^i. A nan¬ 
odiamond hosting multiple NV centres is attached to the 
tip of an AFM. A microscope objective lens placed below 


the sample, which sits on a thin transparent substrate, 
enables us to excite the NV centres and collect their fluo¬ 
rescence. In addition, a microwave antenna is fabricated 
on top of the substrate to allow the NV spin resonances 
to be driven. The quantum probes are type-lb nanodia¬ 
monds ~ 100 nm in diameter, containing N ~ 100 NV 
centres per particle. This leads to a collected fluores¬ 
cence signal of up to 20 x 10 6 photons per second (see 
SI) and a gain in signal-to-noise ratio of y^N « 10 over a 
single NV centre probe. By monitoring the fluorescence 
intensity while sweeping the microwave frequency, a spin 
resonance spectrum of the NV ensemble is obtained (Fig. 
IT)- It generally exhibits two lines at frequencies and 
Z 7 + . By virtue of Zeeman, Stark and thermal expansion 
effects, v_ and are well characterised functions of the 
local magnetic field B [2j [7], electric field E [ 8 ] and tem¬ 
perature T mu* While spin resonance spectroscopy is 
sensitive only to the static component - or sub-kHz vari¬ 
ations - of B, E and T, methods have been developed 
to probe fast fluctuations of these quantities 0 mum. 
These methods rely on quantum techniques to control 
and measure the spin dynamics of the NV centres, which 
are strongly dependent on fluctuations in the local envi¬ 
ronment. Fig. []J illustrates two examples of such mea¬ 
surements, performed on a nanodiamond attached to an 
AFM tip, termed spin relaxation and spin decoherence 
measurements, respectively. The spin relaxation time T\ 
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(left-hand graph in Fig. [1J) is measured by applying 
laser pulses to initialise and subsequently read out the 
NV spin state after a wait time r. The relaxation rate 
Ti si l/Ti can be written as the sum Ti = Ti^t + Fi ?ext 
where Ti^ nt is the intrinsic relaxation rate and Ti^xt is 
the relaxation caused by the fluctuations of the external 
magnetic field occurring at the spin resonance frequen¬ 
cies v± « 3 GHz 11 . [291433] . The spin coherence time 
X 2 (right-hand graph in Fig. El) is measured by apply¬ 
ing MW pulses according to a spin echo sequence. It 
is sensitive to slower magnetic field fluctuations, in the 
kHz-MHz range mmm- In the following, we will use 
spin resonance spectroscopy to image static magnetic and 
temperature fields, but also spin relaxometry to image 
fast magnetic fluctuations. 

Owing to the gain in the collected fluorescence signal 
of over two orders of magnitude compared to the case 
for a single NV centre probe, it is possible to conduct 
imaging experiments with greatly improved acquisition 
speed. Of particular interest is the use of NV quantum 
sensing to characterise the magnetic properties of individ¬ 
ual magnetic nano-objects. As a demonstration experi¬ 
ment, we image and investigate the dynamics of small 
aggregates of maghemite ( 7 ^ 203 ) nanoparticles. Mag¬ 
netic nanoparticles have been exploited in a broad range 
of applications including magnetic seals and inks, mag¬ 
netic recording media, catalysts, and ferrofluids, as well 
as in contrast agents for magnetic resonance imaging and 
therapeutic agents for cancer treatment [34]. Maghemite 
nanoparticles in particular exhibit a rich variety of tun¬ 
able magnetic phenomena, ranging from superparamag¬ 
netism below a critical diameter (^ 10 nm at room tem¬ 
perature) to ferromagnetism for bigger particles and even 
superferromagnetism for some ordered assemblies of par¬ 
ticles [34] . However, it remains a challenge to study the 
magnetic properties of individual - or small clusters of - 
such nanoparticles, especially in their superparamagnetic 
phase [351 [36 b 

Small aggregates of maghemite nanoparticles of diam¬ 
eter 5 — 25 nm were prepared on a glass cover slip (see 
SI). Fig. [ 2^1 shows an AFM image of the sample acquired 
with a nanodiamond on the tip, revealing aggregates of 
various sizes ranging from 50 nm to more than 500 nm, 
corresponding to 10 — 100 particles per aggregate. In 
addition to the topographic image, the fluorescence sig¬ 
nal from the NV ensemble is monitored under a given 
sequence of laser and microwave pulses, and analysed to 
form the corresponding NV sensor image. For instance, 
Fig. |2 Ip shows the fluorescence intensity under continu¬ 
ous laser excitation, without microwave excitation. This 
provides a near-field optical image of the sample, which 
correlates well with the topography image. 

The magnetisation of each magnetic particle can be 
decomposed as M (£) = Mo + $M(£), where Mo is the 
static component and the fluctuating one. Like¬ 

wise, the magnetic field produced by M (£) at the nan¬ 
odiamond location is written as B (£) = Bq + 

To gain information about the magnetic state and dy¬ 


namics of the particles, we applied two different sens¬ 
ing schemes, which rely on spin resonance spectroscopy 
to infer Bq and spin relaxation to infer In Fig. 

Et two MW frequencies v\ and v 2 are consecutively ap¬ 
plied and the normalised difference, the spin resonance 
shift, £>(^ 1 ,^ 2 ) = [^(^ 2 ) ~ E{y\)\lT{y 2 ) is plotted, with 
T(v\^) being the NV fluorescence intensity under MW 
frequency zq ? 2 - Due to the Zeeman effect that shifts 
the spin resonance lines in the presence of a magnetic 
field, this method provides a map of the static compo¬ 
nent Bq of the magnetic field produced by the sample 
[2|. For a single NV centre, the resonance frequencies are 
v± = D dz ybE 2 + (7 e F>o 5 ||/27r) 2 where D and E are the 
zero-field splitting parameters, 7 e /27r = 28 GHz/T is the 
electron gyromagnetic ratio and H 0? y is the magnetic field 
projection along the NV centre’s symmetry axis. For an 
ensemble of NV centres, the spectral features are broad¬ 
ened because the probe experiences four different mag¬ 
netic field projections at the same time corresponding to 
the four possible crystallographic axes, and due to vari¬ 
ations in the zero-field splitting parameters D and E [4|. 
It is seen in Fig. S' that the aggregates produce a mea¬ 
surable static magnetic field component. Fig. shows 
spin resonance spectra acquired with the tip positioned 
at three different locations on the sample, indicated by 
crosses in Fig. EK The spectra are fitted to a model that 
assumes a normal distribution of magnetic field ampli¬ 
tudes B 0j || with a mean value B 0 ,|| an( 4 standard deviation 
c tb 0 || (see SI). The width <tb 0 y accounts for the averag¬ 
ing over the various orientations and positions of the NV 
centres contained in the nanodiamond, as well as over 
the vertical oscillation of the AFM tip. We found these 
values to reach up to H 0 , y = 1.2 mT and gb q y = 0.5 mT. 
This indicates the presence of maghemite nanoparticles 
in a ferromagnetic state, as expected for particles bigger 
than « 10 nm. 

To probe the dynamics of the magnetisation, we mea¬ 
sured the spin relaxation time of the NV ensemble. NV 
relaxometry has been recently applied to the detection of 
paramagnetic molecules [25[ [291133b and superparamag¬ 
netic nanoparticles [H [37]. Fig. [2b presents a spin re¬ 
laxation image formed by recording the fluorescence sig¬ 
nal J-(t) from the nanodiamond after two different wait 
times r\ and 72 , and plotting the normalised difference 
P(ri,r 2 ) = — T(t 2 )\!T(t\). A strong decrease 

in relaxation time is observed for most aggregates. This 
is confirmed by the full relaxation curves measured at 
different locations (Fig. [2]f). The relaxation time is de¬ 
creased from Ti = 39.7 =b 9.2 /is on the bare substrate 
down to T\ = 6.5 db 0.6 //s for the biggest aggregate in 
the investigated area. As the T\ relaxation time of the 
NV centres is governed by the local field fluctuations at 
the resonance frequencies of v± « 3 GHz, our obser¬ 
vations suggest that the maghemite aggregates exhibit 
strong magnetisation fluctuations in the GHz range. As¬ 
suming that the fluctuations of can be described 

by an Ornstein-Uhlenbeck process with a relaxation time 
r m specific to each particle, the induced relaxation rate 
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FIG. 2 . Magnetic imaging of small aggregates of maghemite nanoparticles, (a) AFM image of aggregates of maghemite ( 7 - 
Fe 20 s) nanoparticles on a glass cover slip, (b) NV fluorescence image, (c) NV spin resonance shift image, corresponding to 
the differential signal X>(zq, v 2 ) (see definition in text). The frequencies v\ — 2850 MHz and — 2865 MHz are chosen to give 
a negative (blue) signal in the absence of magnetic held, and a vanishing or positive (red) signal in the presence of a sizeable 
magnetic held, as illustrated in (d). (d) Optically detected spin resonance spectra recorded at three different locations as 
indicated by the crosses in (a) with corresponding colors. The solid lines are the ht to the model described in the text, (e) NV 
spin relaxation image, corresponding to the differential signal T ) (ti,T 2 ) (see dehnition in text). The wait times are n = 1 /is 
and 72 = 20 /xs, as illustrated in (f). (f) Spin relaxation curves recorded at the same locations as in (d). The solid lines 
are the ht to a stretched exponential function, from which we infer the spin relaxation time T±. (g) Schematic illustration of 
superparamagnetism. The magnetic anisotropy energy of a particle is plotted as a function of the angle 0 between the magnetic 
moment and the easy (anisotropy) axis of magnetisation. Thermal fluctuations can cause fast reversals of the magnetic moment 
by overcoming the energy barrier. 


is ri ;ex t ~ (5B(t) 2 ), where (5B(t) 2 ) is the variance 

of the magnetic held amplitude averaged over the NV 
ensemble produced by the particles satisfying the reso¬ 
nance condition r m « (see SI). The measured NV 
relaxation time o f T\ ~ 6 .5 /is therefore translates into a 
held fluctuation (5B(t) 2 ) « (270 /iT) 2 . 

Such fluctuations are consistent with a superparam- 
agnetic behaviour (Fig. [ 2 ^). Indeed, for particles be¬ 
low a certain size, thermal fluctuations can overcome 
the energy barrier due to magnetic anisotropy, result¬ 
ing in spontaneous reversals of the magnetisation direc¬ 
tion. For uniaxial anisotropy, the superparamagnetic re¬ 
laxation time is given by r m = To exp ) ’ w ^ ere K 

is the anisotropy constant, V the particle’s volume, k,B 
the Boltzmann constant and T the temperature. For 
maghemite particles, the prefactor To lies in the range 
HT 11 —1(T 9 s while K « 2-7xl0 4 J/m 3 j39}SI]. For the 
smallest particles with a diameter of just 5 nm, this leads 
to t m = 10“ 11 — 10 -9 s at room temperature, which is 


well within the sensitivity window of NV relaxometry as 
2 ^ 7 } ~ 6 x 10 -11 s. On the other hand, particles as large 
as 25 nm in diameter yield r m > 1 s, which correspond to 
the quasi-static, ferromagnetic state. Our images there¬ 
fore suggest that most aggregates are comprised of both 
small superparamagnetic particles and larger ferromag¬ 
netic particles. Furthermore, from the relative strength 
of the static and fluctuating magnetic signals, we deduce 
that the superparamagnetic particles constitute about 
10% (by volume) of the biggest aggregates in Fig. Ht 
(see SI). We note that this experiment closely follows 
the initial idea in Ref. m to use a scanning NV probe 
to obtain simultaneous static and dynamic information 
(through Ti) about a sample. Measuring the T 2 coher¬ 
ence time would provide further access to low-frequency 
field fluctuations. However, no significant variation in T 2 
time was observed in our experiments, which is due to the 
intrinsic T 2 time being much shorter than the intrinsic T\ 
time [29] . 

An important advantage of using a multiple-NV sen- 
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FIG. 3. Imaging a few-particle maghemite aggregate, (a) 
AFM image of an aggregate made of a few maghemite 
nanoparticles. The area corresponds to the dashed square 
in Fig. [2^. (b) NV fluorescence image, (c) NV spin reso¬ 
nance shift image, (d) NV spin relaxation image, (c) and (d) 
are recorded with the same parameters as in Figs. [2J: and [2^. 
(e) Line cuts extracted from (a), (c) and (d), taken along the 
dashed line shown in each corresponding image. The averag¬ 
ing width is 2 pixels. 


sor is the decreased acquisition time due to increased col¬ 
lected photon signal. The magnetic images shown in Fig. 
Ef and [2^ have been acquired with dwell times of 0.4 and 
1 s per pixel, respectively, corresponding to a total frame 
time of 60 and 150 minutes. This allows us to acquire the 
data in a single pass as in conventional AFM imaging, 
as opposed to the significantly longer acquisition times 
typically needed in single NV imaging experiments [25 1 
128] , which require complex multi-pass procedures and re¬ 
alignment algorithms. We stress that although our NV 
ensemble has inferior features compared with a single 
NV centre in a high-purity bulk diamond (e.g., broad¬ 
ened resonance lines, shorter intrinsic relaxation time) 
which reduces the sensitivity to small changes - i.e., a 
resonance shift much smaller than the line width or a 


relative T\ change much smaller than unity -, this does 
not affect the present experiments where large changes 
are observed. Thus, our images fully benefit from the 
y/N increase in signal-to-noise ratio. Apart from the in¬ 
creased acquisition speed, a further important feature of 
our microscope is its spatial resolution. By imaging a 
cluster of a few (~ 10) maghemite nanoparticles, whose 
topography is shown in Fig. EK the spatial resolution 
of the instrument can be estimated. Figs. HP, hk and 
P show the near-field optical, spin resonance shift and 
spin relaxation images, respectively, obtained with the 
same nanodiamond probe as in Fig. [2| The total acqui¬ 
sition time was 8 minutes for Figs. [3 and [3]:, and 42 
minutes for Fig. P- Again, a fast magnetic fluctuation 
component is observed (Fig. [3]i), while the static field 
seems comparatively weaker (Fig. io- Detailed anal¬ 
ysis of the signal strength confirms that this aggregate 
is mainly composed of superparamagnetic particles 
SI). A line cut across the aggregate is shown in Fig. 

It reveals a feature in the relaxation data as small as 60 
nm (full width at half maximum). This spatial resolu¬ 
tion is limited by the effective probe volume, which is 
given by the size of the nanodiamond (~ 100 nm). We 
note that the magnetic signal is largest on the sides of 
the aggregate. This is due to the fact that the AFM 
tip has a small vertical oscillation, therefore the effective 
probe-sample distance is larger when the tip is above the 
particles - hence a smaller mean magnetic field - than 
when it is positioned at the side of the particle. The 
asymmetry is attributed to the non-spherical shape of 
the nanodiamond. 

Quantum sensing based on NV centres in diamond is 
not limited to the measurement of magnetic fields. It 
has recently been shown that the NV spin can be used 
as a sensitive thermometer by relying on the tempera¬ 
ture dependence of the crystal field splitting parameter 
D [TUI I42H45] . Here we extend this work and demon¬ 
strate the concept of nanoscale temperature mapping. 
As a test sample, we use gold nanoparticles of diameter 
30 ± 10 nm, which exhibit a surface plasmon resonance 
at around 530 nm (see SI). The particles can therefore 
be efficiently heated using the same laser as that used to 
excite the NV centres. Under these conditions, the gold 
nanoparticles act as nanoscale heat sources. A sample 
was prepared by spin coating gold nanoparticles onto a 
glass cover slip (see SI). Because the glass substrate has a 
thermal conductivity much greater than that of air, any 
heat generated by illumination of the gold particle in air 
would be expected to spread preferentially towards the 
substrate, resulting in a negligible change in the air tem¬ 
perature above the particle. Hence, in order to probe this 
photo-induced heating, the nanodiamond probe and the 
sample were immersed in water, as schematically shown 
in Fig. EK Under these conditions, the heat spreads 
quasi-isotropically around the particle, enabling the nan¬ 
odiamond sensor to report the local temperature change. 

In order to locate a gold particle and scan the region 
around it with the nanodiamond probe tip, we relied on 
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FIG. 4. Thermal imaging of a photo-heated gold nanoparticle, (a) Schematic of the experiment. To ensure good thermal 
conductivity between probe and sample, both AFM tip and sample are immersed in water. (b,c) Fluorescence image (b) and 
temperature map (c) obtained simultaneously by scanning a 40-nm gold particle relative to the nanodiamond probe and its 
excitation laser. The integration time is 1.5 s per pixel, hence 50 minutes for the whole frame, (d) Optically detected spin 
resonance spectra corresponding to three different pixels of the scan, located as indicated by the crosses in (b) with matching 
colours. The solid lines are the fit to a sum of two Lorentzian functions. The vertical dashed lines indicate the splitting 
parameter D deduced from the fit, from which we infer the temperature increase AT (see text). The spectra are offset for 
clarity, (e) Line cuts extracted from (b) and (c) taken along the dashed lines shown in (b) with matching colours. The averaging 
width is 4 pixels. The black solid line is a fit to Eq. 0 with the mean diamond-gold particles distance h g d and maximum gold 
particle temperature AT™ ax as free parameters (see text). 


the collected fluorescence rather than on the topography. 
The gold particles produce a weak but measurable flu¬ 
orescence signal, which is maximised when the nanodi¬ 
amond is directly above the particle as the laser beam 
is always focused on the nanodiamond (see Fig. IT)- 
Fig. (4 Jd shows a fluorescence image obtained by scanning 
the nanodiamond around a 40-nm-diameter gold parti¬ 
cle (see SI). The fluorescent spot corresponds to the gold 
particle and indicates a significant absorption of the laser 
light, which is expected to turn primarily into heat. To 
form a temperature map, a full spin resonance spectrum 
is recorded at each pixel of the scan. The spectrum is 
acquired under continuous laser excitation, with an inci¬ 
dent power of 250 fi W, in order to maximise the heat¬ 
ing. From the fits to the spectroscopy data, we then 
infer the splitting parameter D and deduce the tempera¬ 
ture increase AT = T — T 0 , where T 0 is the temperature 
measured in the absence of laser heating (see SI). The 
resulting temperature map is shown in Fig. &• it re¬ 
veals a region of increased temperature up to « 40 K, 
which is spatially correlated with the fluorescence spot 
(Fig. 03). Individual spectra corresponding to a few se¬ 
lected pixels are shown in Fig. Hi. A negative shift of 
the resonance frequencies is clearly observed at the hot 
spot. A line cut across the particle extracted from the 
temperature map is shown in Fig. [4^ (red data), indi¬ 
cating a peak of full width at half maximum ~ 180 nm. 
Also shown in Fig. 03 is a line cut taken far from the par¬ 


ticle (blue data). A small amount of background heating 
of ATbkg ~ 9 K is evident, suggesting that the nanodi¬ 
amond is heated by the laser even in the absence of the 
gold particle, most likely through heating of the silicon 
AFM tip. We stress that the increase in fluorescence 
signal caused by the use of an NV ensemble is crucial to 
this particular experiment. Indeed, the signal from a sin¬ 
gle NV centre would be overwhelmed by the background 
fluorescence produced by the gold particle, making the 
measurement highly challenging using the conventional 
single NV approach. 

To gain more insight into the heating of the gold par¬ 
ticle, we derived an analytical expression for the nan¬ 
odiamond temperature A Td as a function of the relative 
position x g of the gold particle (see SI). The tempera¬ 
ture increase of the gold particle is denoted A T g and is 
assumed to be proportional to the laser intensity at po¬ 
sition x g , while the temperature distribution outside the 
gold particle is obtained from the thermal diffusion equa¬ 
tion in the steady-state regime, assuming for simplicity 
a spherical particle in a homogeneous medium [46] . This 
yields 


^Td(Xg) 


AT bkg + AT™ ax exp 



Rg 



(!) 

where AT™ ax is the temperature of the gold particle 
when it is in the centre of the laser beam (i.e., x g = 0), 
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cjl y/S In 2 = 300 nm is the laser beam width, R g = 20 
nm is the gold particle radius, and h g d is the centre-to- 
centre distance between the gold particle and the nan¬ 
odiamond probe at x g = 0 (see Fig. |4^i) . Fitting Eq. (fl]) 
to the data gives h g d = 70 ± 20 nm, consistent with the 
known sizes of the gold and diamond particles (radius of 
20 and 50 nm, respectively), and AT™ ax = 110 ± 30 K 
(Fig. [4^). This temperature increase can be compared 
to the theoretical value AT™ ax = P a bs/ where 
P abs is the absorbed power, which applies to the case of 
an isolated particle in a homogeneous medium of thermal 
conductivity k. Using the bulk thermal conductivity of 
glass k = 0.8 W-m _ 1 -K _1 and a power P a b s = 9 /YW 
calculated from Mie scattering theory m , one predicts 
AT^nax = 40 K. This is smaller than the value of 110 
K inferred from the experiment, suggesting there may 
be enhanced absorption by the gold nanoparticle in the 
experiment due to, e.g., a non-spherical particle shape 
or tip-induced field concentration effects. Nevertheless, 
these data illustrate that our technique can be used to 
investigate and measure nanoscale heating and thermal 
gradients induced by light, chemical reactions or other 
energy dissipation sources. In addition, Eq. 0 shows 
that the linewidth of 180 nm is a convolution of both the 
temperature profile above the gold particle at a given 
position of the laser beam and the profile of the stimu¬ 
lating laser beam itself. This width is therefore related 
to the measurement conditions rather than to the intrin¬ 
sic spatial resolution, which is limited only by the probe 
volume of ~ 100 nm corresponding to the nanodiamond 
size. We note that direct mapping of the temperature 
profile during constant heating could be achieved by em¬ 
ploying an additional heating laser, thereby providing a 
way to investigate photo-induced heating with a spatial 
resolution drastically better than the current, diffraction- 
limited techniques [48] [49]. 

In this work we have demonstrated a new approach 
to scanning quantum probe microscopy that relies on a 
spin ensemble in a nanodiamond. In addition to con¬ 
ventional AFM topography, our microscope is capable of 
providing spin resonance and spin relaxation images with 
typical acquisition times of an hour. Our instrument is 
capable of imaging both the static and fluctuating mag¬ 
netic fields generated by small aggregates of maghemite 
nanoparticles, enabling detection of both ferromagnetic 
and superparamagnetic particles. The nano-spin ensem¬ 


ble is also shown to be a sensitive nanoscale thermome¬ 
ter providing detailed temperature maps of photo-heated 
gold particles. 

The present approach and experimental apparatus is 
particularly suited to the study of biological samples 
for several reasons. First, scanning nanodiamond mi¬ 
croscopy is compatible with operation in water, which 
paves the way to imaging experiments on living samples 
in their native environment. Second, the use of multi¬ 
ple NV centres provides a significant gain in acquisition 
speed, which is a clear advantage for fragile samples that 
cannot tolerate long acquisition times. In the experi¬ 
ments reported here the nanodiamond contained ^100 
NV centres, but even higher concentration of NV centres 
could be used to further improve the scanning speed [50] . 
Lastly, our instrument provides topographic images cor¬ 
related to the quantum sensor images, with a combined 
spatial lateral resolution below 100 nm. This is a crucial 
requirement for experiments performed on samples with 
a non-trivial shape. We envisage using this technology to 
investigate spin signals in cell membranes stemming, e.g., 
from spin labels [3T or ion channels EU, which could give 
crucial insight into cell membrane functions and dynam¬ 
ics. 

On the other hand, scanning quantum probe mi¬ 
croscopy also opens up opportunities for the study of 
solid-state samples. In particular, the temperature map¬ 
ping capability reported in this paper provides a new 
platform for studying thermal transfer at the nanoscale. 
This is particularly relevant in the race towards ever 
smaller micro/nano-electronic devices, where the precise 
characterisation of heat transfer remains a fundamental 
bottleneck [52] . 
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SUPPORTING INFORMATION 


I. EXPERIMENTAL SETUP 

The experimental apparatus is based on an AFM (Asylum Research MFP-3D-BIO) mounted on a customised 
inverted optical microscope (Olympus IX-71). The microscope objective lens (Olympus PlanApoN 60x, NA = 1.42) is 
mounted on a 3-axis XYZ scanning stage (PI P545-3R7) to allow fast laser scanning. The filter cube (Olympus BX 2 ) 
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FIG. 5. (a,b) Confocal images of the AFM tip in the vertical YZ plane (a) and horizontal XY plane (b), the XYZ reference 
frame being defined in main-text Fig. 1. The dashed (dotted) lines depict the AFM tip (the sample’s top surface). The bright 
spot at the apex of the tip corresponds to the grafted nanodiamond. In (a), the left-hand (right-hand) image shows the tip 
retracted (engaged) relative to the sample’s surface, (c) Number of photons collected per second P(Pl) from a nanodiamond 
on tip as a function of the incident laser power Pl- The solid line is the fit to a saturation law P(Pl) — 1+ p°l/p L ? yielding 

the fit parameters JTo = 33.3 dfc 1.4 X 10 6 counts/s and P sa t = 275 =b 26 fi W. 


loaded in the microscope turret is equipped with a dichroic beam splitter (Chroma ZT532rdc) and an emission filter 
(Chroma NC359095 ET690/120x). The excitation laser (Laser Quantum Gem 532) is modulated using an acousto¬ 
optic modulator (AA Opto-Electronic MT-200-VIS) in a double pass configuration. The modulated laser beam feeds 
the back port of the Olympus IX-71 microscope through a single-mode fibre. The filtered fluorescence light exits the 
microscope through the side port and is coupled into a multi-mode fibre connected to a single photon counting module 
(Excelitas SPCM-AQRH-14-FC). The photon signal is analysed by multi-function devices (NI-DAQ PCIe-6323 and 
PCIe-6321) for simple counting and a time digitizer (FastComTec MCS6A-2T4) for time-correlated measurements. 

The sample is mounted onto a 2-axis XY scanning stage while the AFM cantilever sits on a Z piezo-electric ac¬ 
tuator for distance control. The NI-DAQ devices are also used to generate analog voltages that control the XYZ 
objective stage and the XY AFM stage, as well as to read out the Z position of the AFM cantilever, which is directly 
regulated by the AFM controller. The sample is equipped with a microwave antenna fabricated by optical lithography 
and metal deposition. It is connected to the output of a microwave generator (Rohde & Schwarz SMB 100A) after 
modulation (Mini-Circuits ZASWA-2-50DR+) and amplification (Mini-Circuits ZHL-16W-43+). The laser and mi¬ 
crowave modulations are controlled by a programmable pulse generator (SpinCore PulseBlasterESR-PRO 500 MHz). 
All instruments and devices are controlled with a single purpose-built Labview program. 


II. PROBE PREPARATION AND CHARACTERISATION 

The nanodiamonds used in this work (rFND-100) were provided by the BioDiamond team of Academia Sinica 
(Taipei, Taiwan). They were spin cast on a glass coverslip and characterised both optically and by AFM. The nanodi¬ 
amonds selected to be attached to the AFM tip were ^100 nm in size and contained between 20 and 100 NV centres 
based on the observed fluorescence intensity. The AFM probes used here are either medium-soft silicon cantilevers 
with a nominal tip radius of 9 nm (Olympus AC240TS-R3) for the magnetic imaging experiments or stiff silicon 
cantilevers with a 2-/rm-diameter plateau tip (Nanosensors PL2-NCHR-10) for the thermal imaging experiments. 

To transfer a nanodiamond to an AFM tip, we must first align the laser beam onto the AFM tip. A coarse alignment 
is achieved using a CCD camera placed on top of the AFM head, providing a precision of ±2 /i m. The alignment is 
then fine tuned by correlating the AFM and confocal images obtained when scanning the sample. Once a nanodiamond 
suitable for sensing - i.e. with a good spin contrast - is found, an AFM scan centred on this nanodiamond is performed 
repeatedly. By decreasing the oscillation amplitude of the tip in AC (tapping) mode or by switching to contact mode, 
the effective interaction strength between the tip and the nanodiamond is increased, enabling the tip to push the 
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(a) Nanodiamond probe 'ND1' used in Figs. 2 and 3 



MW frequency, v (MHz) 

(b) Nanodiamond probe 'ND2' used in Fig. 4 




MW frequency, v (MHz) 





Echo time, 2r (ps) 


FIG. 6. Characteristic curves measured for ND1 (a) and ND2 (b). The graphs in the left, middle and right columns correspond 
to optically detected spin resonance spectra, spin relaxation curves and spin echo decoherence curves, respectively. The markers 
are the data point while the solid lines are fits as described in the text. 


nanodiamond on the substrate and eventually attach it [T9] . The nanodiamond grafting is evidenced by a strong 
increase in the fluorescence signal from the tip - the laser being aligned on the tip at all times - and an increase in 
the tip height. To test and strengthen the bonding, we then perform long, low-speed scans in AC mode. 

Figs. and show confocal images of an AFM tip after grafting a nanodiamond onto it. The nanodiamond 
appears as a bright fluorescence spot at the apex of the tip. Fig. [5j^ shows a typical saturation curve from this spot, 
that is, the collected fluorescence signal as a function of the incident laser power. In this example the signal reaches 
up to 22 x 10 6 counts/s close to saturation, with a laser power of 600 /iW. Once a nanodiamond is attached to the 
AFM tip, its performances as a quantum sensor can be characterised. Fig. [6] shows the relevant characteristic curves 
for the nanodiamond probes used in this work, the top row corresponding to the probe ‘NDl’ used in main-text Figs. 
2 and 3 for the magnetic imaging experiments (Fig. it), while the bottom row corresponds to the probe ‘ND2’ used 
in main-text Fig. 4 for the thermal imaging experiments (Fig. if). The left-hand column in Fig. [6] corresponds to 
the spin resonance spectra T[v\ i.e. the fluorescence intensity as a function of the MW frequency. They are fitted to 
a sum of two Lorentzian lines 


T{y) = To 


1 - 


C- 


1 + 


Au- 


2 


i + 


C± 

A^_)_ 


( 2 ) 


where To is the off-resonance fluorescence intensity, and z/±, C± and Av± are the central frequency, relative contrast 
and full width at half maximum, respectively, of each of the two spin resonances denoted by the subscript ±. In the 
absence of external magnetic field, the resonance frequencies are v± = D ± E, D and E being the so-called zero-field 
splitting parameters of the NV centre. Note that since the nanodiamonds contain an ensemble of NV centres, the 
spectra indicate the mean characteristics of this ensemble. Here we find D = 2868.9 ±0.1 MHz and E = 7.5 ±0.1 
MHz for ND1 and D = 2868.0 ±0.2 MHz and E = 5.3 ± 0.2 MHz for ND2. Under typical experimental conditions, 
the contrast C± is about 8% and the linewidth is A u± ~ 10 MHz. 

The middle column in Fig. [6] is the spin relaxation curve T(r) obtained by measuring the fluorescence intensity 
following an initialisation laser pulse and a wait time r. The resulting curves are fitted to a stretched exponential 
function 


Hr) = 


1 -C + Ce ±) 


(3) 


where To is the fluorescence intensity in the initialised spin state, C is the relative contrast, T\ is the spin relaxation 
time of the NV ensemble and n is an exponent that accounts for the inhomogeneity of the spin ensemble. In our 
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nanodiamonds we have typically C ~ 15%, n « 0.5 — 1 and T\ « 40 — 400 /is. The Ti values of ND 1 and ND 2 are 
indicated in Fig. [ 6 | 

Finally, The right-hand column in Fig. [ 6 ] is the decoherence curve obtained under a spin echo sequence ^ — r — 
7r — r — The resulting curves are fitted to a single exponential function J r (2r) = ^[1 — C + C exp(— 2t/T 2 )], with 
X 2 being the spin coherence time. In our nanodiamonds, X 2 is typically comprised between 500 ns and 1.5 /is in zero 
external magnetic field. 


III. SUPPORTING INFORMATION FOR THE MAGNETIC IMAGING EXPERIMENTS 


A. Sample preparation and characterisation 


The magnetic nanoparticles used for the experiments reported in main-text Figs. 2 and 3 were purchased from 
Sigma-Aldrich (reference 544884). They are received in the form of a nanopowder of Fe 2 C >3 primarily in the 7 
crystalline phase, or so-called maghemite. We diluted the powder in Milli-Q water to obtain a solution of concentration 
of 0.1 g/L. After plasma cleaning a glass cover slip for 10 minutes, 50 /iL of the maghemite solution is dropped on the 
cover slip and immediately spinned using a spin caoter with an acceleration of 200 rpm/s and a speed of 3000 rpm for 
60 s. These parameters result in a dispersion of particles on the substrate that are mostly isolated from each other, 
with only little aggregation, as shown in the AFM image in Fig. Ob The height of a particle in the AFM image gives 
an estimate of its mean diameter (Figs. [7 ]d and Ob- We find that the diameter of the particles in this sample ranges 
from 3 nm to 25 nm. This is in good agreement with the surface area of 50 — 245 m 2 /g specified by the manufacturer, 
which corresponds to a diameter of 5 — 25 nm. To obtain aggregates of the particles, we repeat the above procedure 
but with a spinning speed decreased to 1500 rpm, thus increasing the probability for aggregation. An AFM image of 
the sample after this step is shown in main-text Fig. 2a. 

We also characterised the magnetic properties of a macroscopic sample of the same maghemite nanoparticles. We 
placed ~ 1.6 mg of the as-received powder in a SQUID magnetometer and measured the average magnetisation M 
as a function of magnetic field H and temperature T. Fig. [ 8^1 shows the M(T ) curve under an applied field of 
/io H = 1 T, which should saturate all particles regardless of their size. The small decrease of M at 300 K relative to 
the low temperature value indicates that the Curie temperature is well above 300 K. Fig. [ 8 Jd shows the M(H) curve 
at T = 300 K, obtained by ramping down the field from 5 T to 0. The remnant magnetisation at H = 0 indicates 
the presence of ferromagnetic-like particles. For an ensemble of ferromagnetic particles with randomly oriented easy 
axes, one would expect a ratio m(h-¥oo) ~ contrary, if the particles were purely superparamagnetic, the 

magnetisation would average to zero over the measurement time, leading to ~ 0- The P resen t data (Fig. 8p) 

suggests that the powder comprises both ferromagnetic-like (in fact, ferrimagnetic) and superparamagnetic particles, 


with approximately equal volumes since 



FIG. 7. (a) AFM image of a sample of non-aggregated maghemite nanoparticles, (b) Close-up view of a single particle, (c) 
Line cut taken along the red dashed line in (b). 
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T (K) 



FIG. 8. (a,b) Characteristic curves of a macroscopic sample (1.6 mg) of the maghemite particles measured using a SQUID 
magnetometer, (a) Magnetisation versus temperature, under an applied field of 1 T. (b) Magnetisation versus applied field at 
300 K. The field is ramped down from 5 T to 0. 


B. Data acquisition methods 


The images reported in main-text Figs. 2 and 3 were obtained by scanning the sample stage while operating the 
AFM tip in AC (tapping) mode. The amplitude of the vertical oscillation was estimated from force curves to be ~ 60 
nm in scanning conditions. The NV fluorescence signal is integrated only once for a given time at each pixel, after 
which the sample is moved to the next pixel. After each line of the scan, the XYZ position of the objective lens is 
adjusted in order to maximise the NV fluorescence signal. For the spin resonance images (main-text Figs. 2c and 
3c), the microwave is switched on for 150 ns, corresponding to a tt flip rotation of the NV spins, followed by a 300-ns 
laser pulse for readout and a l-/is wait time. This sequence is repeated continuously during the pixel time, with the 
microwave frequency set to zq for half the pixel time and zq for the second half. The integrated fluorescence signals 
J*(z/i) and T(y 2 ) a re recorded using the NI-DAQ device. For the spin relaxation images (main-text Figs. 2e and 3d), 
a sequence of three 3-/is-long laser pulses with wait times t\ and r 2 is repeated continuously during the pixel time. 
The fluorescence signal is acquired using the FastComTec time digitizer and the first 300 ns of each laser pulse are 
integrated to give the values T[t\) and T{j 2 ). In these experiments, the incident laser power was Pl ~ 150 /iW. 
All measurements were performed at room temperature under normal atmosphere, with no external magnetic field 
applied. 


C. Analysis of the spin resonance spectra 

Here we describe the model used to analyse the spin resonance spectra recorded on the maghemite aggregates 
(main-text Fig. 2d). The static component of the magnetic field at the location of NV centre i is denoted Bq. This 
NV centre will exhibit spin resonances at frequencies 


4 


= D± 



Pe TDl 

2tt °’H 


(4) 


where B * y = Bq • e\ is the magnetic field projection along the NV centre’s symmetry axis, defined by the unit 
vector e l z . Experimentally, one probes a ensemble of NV centres with different symmetry axes and positions within 
the nanodiamond of diameter « 100 nm. In addition, the AFM tip is operated in tapping mode, which implies the 
existence of a vertical oscillation of ~ 60 nm at a frequency of ~ 60 kHz. These effects all contribute to produce a 
broad distribution of field projections B q y, resulting in a broad distribution of frequencies zdj_. We characterise this 

distribution by a mean value B 0j || and standard deviation gb q y . That is, the probability of having a field projection 
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F> 0 ,|| in the NV ensemble during the measurement is given by 


V{B U ) 



(J3 o,ll 

2 cr 



2 


( 5 ) 


The spectra in main-text Fig. 2d were fitted by a superposition of spectra in the form of Eq. © with frequencies 
given by Eq. Q and with the probability distribution (jHl. The only fit parameters are E> 0 ,|| and ctb q y • All other 
parameters (D, E, C± and A is±) are extracted from a reference spectrum recorded far from any magnetic particle 
(black curve in main-text Fig. 2d). 

It is useful to relate the mean field projection H 0 ,|| to the mean field amplitude Bq. Let us write the field projection 
seen by a given NV centre as B * y = Bo cos 0^, with Oi being the angle between the local field Bo and the NV centre’s 
axis. Assuming that all directions are equiprobable within the Air solid angle, we obtain 


The maximum field measured while scanning the sample of maghemite aggregates corresponds to the spectrum shown 
in main-text Fig. 2d (green curve). The fit to the data yields a mean projection F> 0? y = 1.2 ± 0.1 mT, hence an 
amplitude Bq = 2.4 db 0.2 mT, with a width gb 0 y = 0.5 ± 0.1 mT. 


D. Analysis of the spin relaxation data 

We now describe the relation between the measured spin relaxation time (7\) and the local magnetic field fluctu¬ 
ations. The fluctuating component of the magnetic field at the location of NV centre i is denoted TB 2 (£). Its spin 
relaxation rate reads m 


1 tot 1 int 


( 7 ) 


where 


r; xt = SB i L (t) 2 )f(27rDr m ). (8) 

Here (...) designates the time average, 5B]_(t) 2 = 5B l x {t) 2 +SBy(t) 2 is the transverse component of the field (z being 
defined by the NV centre’s symmetry axis) and r m is the correlation time of The function f{pc) is defined as 


f( x ) 


2x 

1 + X 2 ’ 


(9) 


which is peaked around x = 1 . Since the field 8B\_(t) is produced by an assembly of magnetic nanoparticles with 
different sizes hence different dynamics r m , the relaxation rate Fg Xt will be dominated by the particles for which 
r m ~ 2 ttD ’ m aximising f{2'KDr m ). One can therefore write 

r --^(^i (i)2) ( 10 ) 

where 5B]_(t) now accounts only for the contribution of the particles satisfying r m ~ 2 nD‘ Averaging over all NV 
centres and assuming that all directions of the NV axis are equiprobable, we obtain an effective relaxation rate for 
the nanodiamond 


r 


ext 


2 nD 


(SBm, 


(ii) 


where SB(t) 2 is the squared amplitude of the field seen by a given NV centre and ~ designates the average over the 
ensemble of NV centres within the nanodiamond. We can therefore infer this amplitude from the measured relaxation 
rate according to 


2t rD 


(SBm 


( 12 ) 
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The maximum relaxation rate measured while scanning the sample of maghemite aggregates corresponds to the 
data shown in main-text Fig. 2f (green curve). By comparing with the reference curve (black curve), one deduces 

rext = Ttot-Tint = 129T15 x 10 3 s -1 . This translates into a root-mean-square (rms) amplitude yj(5B(t) 2 ) « 270 fiT- 
rms. 

Finally, we can use the ratio between to the static field amplitude Bq (up to ~ 2.4 mT in main-text Fig. 2c) and the 
fluctuating component ( 5B(t ) 2 ) to estimate the relative populations of ferromagnetic and superparamagnetic particles 
in the sample. Let us express the effective magnetic moment of an aggregate of particles as m(£) = m 0 + Sm(t), 
where mo is associated to the ferromagnetic particles and Sm(t) to the fast-fluctuating superparamagnetic particles. 
The magnetic field produced by the aggregate writes B (t) = Bq + <SB(t). Because of spatial averaging over many NV 
centre locations within the nanodiamond, it is a good approximation to write 


yj ( Sm(t ) 2 ) ^ \J 
m 0 ~Bq 


(13) 


Since the magnetic moment m is proportional to the volume of the particles, the ratio (13) can be 
ratio in volume between superparamagnetic and ferromagnetic particles, which is of the order of 
the data of main-text Fig. 2. This ratio is consistent with the measured particle size distribution 
the estimated threshold size of 10 nm for the ferromagnetic/superparamagnetic transition. 


interpreted as the 
10% according to 
of 5 — 25 nm and 


E. Supplementary data 


In main-text Fig. 3, spin resonance shift and spin relaxation images of a small aggregate of maghemite particles 
are shown (AFM image reproduced in Fig. [9k). These images represent the differential signal £>(^ 1 ,^ 2 ) = [Ffa) — 
F(v\)\ jB( 1 * 2 ) and D(ti, 72 ) = [J~ (ti ) — T (t 2 )\ / T (ti) , respectively. To quantify more precisely the strength of the 
magnetic field produced by the sample, we show here full spin resonance spectra (Fig. §0 and Ti spin relaxation 
curve s (Fig. [9^) recorded at two different locations on the sample. Following the same analysis as in sections III C and 

we find a maximum value of \j (5B(t) 2 ) « 160 /iT-rms for the fluctuating component, while the static component 


HID 


is Bq « 100 fiT. We conclude that the aggregate imaged in main-text Fig. 3 is dominated by superparamagnetic 
particles. 



FIG. 9. (a) AFM image of a small aggregate of maghemite particles, reproduced from main-text Fig. 3a. (b) Optically detected 
spin resonance spectra recorded with the tip at two different positions, indicated by crosses in (a). The solid lines are data fit 
as explained in Section [III C| T he black curve is a reference spectrum recorded far from the aggregate, while the orange curve 
close to the aggregate gives Bo = 0.1 ± 0.1 mT and cfb 0 m = 0.0 ± 0.1 mT. (c) Spin relaxation curves recorded at the same two 
positions. The solid lines are data fit to a stretched exponential function, yielding a reference relaxation time T\ — 51.3± 16.2 fi s 
(black curve) and T\ = 15.8 ± 2.3 /xs (orange curve). 


F. Theoretical estimates of the magnetic field produced by a magnetic particle 

In this section, we compare our measurements of the magnetic field produced by aggregates of maghemite particles 
to simple theoretical estimates. In the experiment, the nanodiamond has a diameter of ~ 100 nm and oscillates 
vertically with an amplitude of ~ 60 nm. This implies that the nanodiamond centre is located on average at a 
distance h ~ 110 nm above the sample’s surface. Since h is significantly larger than the typical diameter d of the 
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maghemite particles (d = 5 — 25 nm), we can treat the particles as point-like magnetic dipoles in order to calculate 
the field probed by the nanodiamond. 

The magnetic moment of a particle of diameter d and magnetisation M has a magnitude m = and is oriented 

along the easy axis of the particle. The amplitude of the magnetic field right above the particle is then given by 

B =^(l + 3c°s 2 9 m ), (14) 


where 0 m is the angle between the moment m and the vertical direction and h is the distance to the particle. Using 
M = 300 kA/m (obtained from Fig. [ 8 ) 3 ), h = 110 nm and d — 20 nm, Eq. (14) predicts that the field amplitude varies 
between 100 fiT for = 7 r /2 and 400 (i T for 0 m = 0. Therefore, an aggregate of several particles is expected to 
produce a field in the mT range, with a static component associated to the ferromagnetic particles and a fluctuating 
one caused by the superparamagnetic particles. This order of magnitude is in agreement with the experiment, where 
field amplitudes up to 2.4 mT (static) and 270 /iT-rms (fluctuating) were measured. 


IV. SUPPORTING INFORMATION FOR THE THERMAL IMAGING EXPERIMENTS 


A. Sample preparation and characterisation 


The gold nanoparticles used for the experiments reported in main-text Fig. 4 were synthesised following the 
protocol described in Ref. [53]. Briefly, we proceed in two steps. First a gold seed solution is prepared, with a 
diameter d = 13 =b 1 nm. To this end 100 mL of a 1.0 mM aqueous HAuCUiSLUO solution was added to 100 mL of 
milli-Q water, which was then boiled. 10 mL of a 38.8 mM aqueous solution of sodium citrate was added and boiled 
for 30 min. Next, a seed-mediated growth approach was applied in order to synthesize gold particles of desired size. A 
volume of 4 mL of a 20 mM aqueous HAuCl 4 : 3 H 2 0 solution and 0.4 mL volume of a 10 mM aqueous AgNCU solution 
were added to 170 mL of milli-Q water. 3 mL of the 13-nm gold seed solultion was added to the solution. 30 mL of 
5.3 mM ascorbic acid was then added slowly (0.6 mL/min) with constant stirring. 


The resulting solution contains gold particles with a mean diameter of « 30 nm, as measured by dynamic light 
scattering. Fig. [lO] shows an extinction spectrum of the solution. It exhibits a surface plasmon resonance at a 
wavelength of ~ 530 nm, which matches well the wavelength of the laser used in the experiment to excite the NV 
centres. This implies that the same laser will be efficiently absorbed by the gold particles, eventually causing significant 
heating. A sample of mostly isolated gold particles was prepared by spin casting a drop of the solution on a glass 
cover slip. An AFM image of the sample is shown in Fig. 11 r. By measuring the height of each particle (Figs. Eyf 
and [Tlfc) , we find the diameter to vary between 20 and 40 nm. 



Wavelength (nm) 

FIG. 10. Extinction spectrum of the solution of gold nanoparticles recorded prior to deposition on the substrate. The surface 
plasmon resonance matches well the wavelength of the excitation laser used in the experiment (532 nm). 
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FIG. 11. (a) AFM image of a sample of gold nanoparticles, (b) Close-up view of a single particle, (c) Line cut taken along the 
red dashed line in (b). 


B. Pre-characterisation of the gold particle investigated 


For this experiment, an AFM cantilever equipped with a plateau-shaped tip was used and the nanodiamond was 
attached to the side wall of the plateau (see Fig. [T2 ^l). As explained below, this enables us to scan the nanodiamond 
above the gold particle without any direct contact between tip and particle, thus preventing unwanted friction effects. 

Fig. [I2| ) shows a confocal image recorded prior to the AFM scans. The position of the nanodiamond as well as the 
plateau tip to which it is attached are indicated. Also shown is the gold particle that was scanned in main-text Fig. 
4. Its fluorescence signal is measurable although significantly weaker than that of the nanodiamond. Figs. [i~2} ' and 
12 i show the AFM and fluorescence images obtained simultaneously by scanning the sample (i.e., the gold particle) 
while focusing the laser on the nanodiamond. Since the AFM image reflects the convolution between the particle 
shape and the tip shape, the particle appears here as a crescent moon due to the plateau being slightly tilted with 
respect to the substrate’s surface. The maximum height observed in the AFM image is indicative of the size of the 
gold particle, which is « 40 nm (Fig. [12^). On the other hand, the fluorescence image exhibits a bright spot stemming 
from the gold particle. Since the laser is focused on the nanodiamond, this spot corresponds to the situation where 
the nanodiamond is exactly above the gold particle. The absence of overlap between the fluorescent spot and the 
topography feature of the gold particle implies that the nanodiamond can be scanned around the gold particle with 
a constant nanodiamond-substrate distance, without direct contact between tip and gold particle. Furthermore, it is 
sufficient to use the fluorescence signal to centre the scan on the gold particle, without resorting to the topography 
information. In main-text Figs. 4b and 4c, we thus scanned the nanodiamond probe around the gold particle while 
recording spin resonance spectra, allowing us to form a temperature map of the photo-heated gold particle. 


C. Acquisition of the temperature map 


The images reported in main-text Figs. 4b and 4d were obtained by scanning the sample stage while operating 
the AFM tip in AC mode. A drop of milli-Q water was added onto the substrate prior to mounting the AFM head. 
We made sure that the AFM cantilever was fully immersed in water by monitoring its resonance frequency, which 
dropped from « 300 kHz in air to ~ 120 kHz in water. The oscillation amplitude was estimated from force curves to 
be ~ 10 nm in scanning conditions. The XYZ position of the objective lens was adjusted after each line of the scan 
in order to maximise the NV fluorescence signal. In main-text Fig. 4d, a spin resonance spectrum composed of 25 
frequency bins was acquired at each pixel, some of the obtained spectra being shown in main-text Fig. 4c. The time 
per bin is 60 ms, resulting in a pixel time of 1.5 s. Both laser and microwave were applied continuously during the 
acquisition, with an incident laser power of Pl « 250 fi W. 

Each individual spectrum was fitted using Eq. ©• The obtained resonance frequencies v± allow us to infer the 
splitting parameter D = u++ ^ . The temperature increase AT plotted in main-text Fig. 4d was calculated according 
to 


A T = 


D-D 0 
-75 kHz/K’ 


(15) 
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Fluorescence 
from the gold 
particle 


FIG. 12. (a) Schematic of the experiment showing the nanodiamond attached to the side wall of a plateau-shaped tip. (b) 
Confocal image obtained by scanning the laser beam relative to sample and AFM tip, with the AFM tip engaged. The red 
circle depicts the 2-/xm-diameter plateau. (c,d) Topography image (c) and fluorescence image (d) obtained simultaneously by 
scanning the sample relative to the AFM tip, with the laser being focused on the nanodiamond. The red circle delimits the 
topographic feature that would results from the plateau being perfectly parallel to the substrate’s surface, (e) Line cut taken 
along the white dashed line in (c), indicating a particle size of ~ 40 nm. 



Laser power (pW) 


FIG. 13. Splitting parameter D of the NV ensemble as a function of the incident laser power Pl . The measurement conditions 
are the same as in main-text Fig. 4c, i.e., the nanodiamond is attached to the plateau tip, the tip is engaged far from any gold 
particle and the laser is focused on the nanodiamond. Shown on the right-hand side is the corresponding temperature increase 
AT deduced from Eq. (15). 


where the numerical factor is the linear conversion factor measured by Acosta et al. near room temperature m- 
The offset Dq is the splitting parameter measured in the absence of laser heating. To obtain Do, we recorded spin 
resonance spectra as a function of the laser power Pl, with the tip engaged on the sample far from any gold particle. 
Fig. [13] shows the parameter D extracted from these spectra as a function of Pl- A linear fit provides us with 
Dq = D(Pl = 0) = 2868.9 db 0.2 MHz and a slope of —2.8 =t 0.4 kHz//iW. This slope indicates the existence of 
a laser-induced heating effect even in the absence of the gold particle. This most likely occurs because some laser 
light is absorbed by the silicon AFM tip, thereby heating the nanodiamond attached to it. With the laser power 
of Pl ~ 250 fi W used in the experiment (main-text Fig. 4c), this intrinsic heating is AT « 9 K. However, this 
temperature increase is significantly smaller than that caused by laser-induced heating of the gold particle, as clearly 
seen in main-text Fig. 4e. 
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D. Derivation of the analytical formula used to fit the measured temperature profile 


In this section we derive the analytical formula given in the main text (Eq. (1)), which represents the nanodiamond 
temperature A Td(x g ) as a function of the position x g of the gold particle relative to the nanodiamond and laser beam, 
as experimentally measured in main-text Fig. 4e. When the gold particle is moved horizontally to a new position x g , 
two effects contribute to change the temperature of the nanodiamond probe: (i) the gold particle receive a different 
amount of heating power because of the intensity profile of the stimulating laser beam, and (ii) for a given heating 
power the temperature at the diamond location is changed because the distance to the gold particle is changed. 
Hereafter these two independent effects are modelled in order to obtain an analytical formula. 

a. Heating power vs. position The transverse profile of the laser beam is modelled by a 2D Gaussian distribution 
with a full width at half maximum (FWHM) of 300 nm determined experimentally through confocal imaging of a 
single NV centre. The power density is denoted as Pl(i*) with r = (x,y) being the transverse position, and reads 

Pi(r,= 2 (16> 


where cr^ a/8 In 2 = 300 nm and Pl = fjP l { y )^ 2r = 250 fi W is the total laser power. We consider a spherical gold 
particle of radius R g . The power absorbed by the particle can be written as 


P&bs{% g) — ^abs^L {Xg) -> 


(17) 


where cr a b s is the effective absorption cross section of the particle and Il{x 9 ) is the laser intensity averaged over the 
cross-sectional area of the particle, that is, 


Il(x 9 ) = —— f [ p L ( r)d 2 r, 

tfgeom J JA(x g ) 


(18) 


with A{x g ) being the geometric cross section of the particle of total area c t { 
20 nm stimulated at a wavelength of 532 nm, Mie theory m predicts a ratio 
absorbed power of P a bs(0) 


geom — 7rR g . 


crabs 

crgeom 


For a gold sphere of radius 
3.0, which yields a maximum 


9.0 fi W. Fig. 14 shows the relative absorbed power P a bs(^p)/-fabs(0) calculated from 


Eq. (18), as well as from the approximation by a Gaussian function 


P&bs (x g ) 
Pahs ( 0 ) 



(19) 


with cjl a/8 In 2 = 300 nm being the FWHM of the laser beam. In the following, we will use the Gaussian approximation 
(19), which is in excellent agreement with the exact integration. Moreover, we will assume that all the absorbed power 
is converted into heat, that is, P au b s {x g ) represents the heating power generated by the gold particle for a given position 
x„. 



Position of the gold particle, (nm) 


FIG. 14. Normalised power absorbed by a gold nanosphere as a function of the relative position of the exciting beam. The 
radius of the sphere is 20 nm and the FWHM of the Gaussian beam is 300 nm. 
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b. Temperature distribution We now analyse the temperature distribution outside the gold particle for a given 
heating power P a bs- Some of the assumptions and approximations made in this paragraph will be justified in the next 
section by comparison with numerical simulations. In the steady-state regime, the thermal diffusion equation is 


V • [«(r)VT(r)] = -p h eat(r), 


( 20 ) 


where «(r) is the position-dependent thermal conductivity and Pheat( r ) is the power density of the heat source. Here, 
the only source of heat is the gold particle with a total power fff pheat( r )d 3 r = P a bs- The particle is modelled as 
a sphere of radius R g in a homogeneous medium of thermal conductivity K e ff- Because the thermal conductivity of 
gold is much larger than that of the surrounding medium, the temperature in the gold particle does not depend on 
the exact distribution Pheat(r) and can be considered as uniform [46] . with a value 


A T n = 


labs 


Outside the gold sphere, the temperature distribution 


AirR g ' 


reads 


. _ -LLn 




( 21 ) 


( 22 ) 


where r > R g is the distance from the gold sphere centre. In the experiment, the gold particle is moved horizontally 
(position x g ) while the nanodiamond probe is standing still at a height h g d , defined as the centre-to-centre distance 
between the gold sphere and the nanodiamond when the two particles are vertically aligned (i.e. x g = 0). We assume 
that the diamond temperature A T d can be approximated by the value calculated at the centre of the diamond volume 
in its absence, that is, 


a T d (x g ) » AT(r = ^/xf+h^ d ) « A T g (x g ) 


Rn 


Using Eqs. (19) and (21), one can write 


JT+K 


(23) 


gd 


Tg( x g) ~ ^g(0 ) ex P ( 2^-2 J ? 


(24) 


where T g ( 0) = P a b s (0)/47rR g K e ff • Inserting Eq. (24) into Eq. (23), one finally obtains 


A T d (x g ) « AT g (0) exp - 


Rn 


2al 


hi 


(25) 


gd 


which is Eq. (1) of the main text. The temperature of the gold particle at x g = 0 (T g (0)) and the height of the diamond 
probe (hgd) are used as fit parameters, while the other parameters are taken to be R g = 20 nm and ctlVS In 2 = 300 
nm. 

c. Comparison with numerical simulations COMSOL Multiphysics (v5.1) was used to simulate the temperature 
distribution of a system comprising a gold nanosphere on silica glass and a diamond nanosphere suspended in water 
(Fig. |15 |l). The simulations were performed on a diamond nanosphere (radius 50 nm) suspended in a hemisphere 
of water (radius 50 /im) 60 nm above a cylinder of silica glass (radius and height 50 yum), and a gold nanosphere 
(radius 20 nm) placed on the surface of the silica glass (Fig. [I5j d). These dimensions were chosen to be representative 
of the experimental situation. The gold nanosphere is heated with a constant power P a bs(^) following Eq. (19). 
As a test value, we use a heating power of P a b s (0) = 36 fi W. A boundary condition of constant temperature at the 
ambient temperature of 296.15 K was set to replicate an infinite medium of water and silica glass. The calculation 
domains were meshed with a physics controlled extremely fine triangular mesh to ensure correct geometric treatment 
of the gold and diamond nanospheres and to increase accuracy of heat transfer. The program then solves the thermal 
diffusion equation (20) in the steady-state regime. We note that the steady state is reached within a /is time scale 
m, which is much faster that our measurement time. 

We first consider the case where the gold nanosphere is aligned with the nanodiamond and the laser beam, i.e. 


x g = 0. The resulting temperature map in the xz plane is shown in Fig. 16 r. Also shown for comparison is the case 


where the diamond is replaced with water (Fig. [l6p). Figs. 16>e show line cuts taken along the vertical direction 
(c), and along the horizontal direction across the gold sphere (d) and across the diamond sphere (e). In the absence 
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40nm 


FIG. 15. (a) Geometry of the domains considered in the thermal simulation, (b) Close-up view showing the dimensions. The 
simulation is run for different positions of the gold nanosphere. For each position x g , a constant heating power Phe&t(x g ) is set. 


of the diamond, the temperature increase inside the gold sphere is nearly uniform and reaches A T g ^ 177 K. Outside 
the gold sphere, the temperature decays quasi isotropically because the thermal conductivity of glass and water are 
similar. The presence of the diamond sphere affects only marginally the general temperature distribution around the 
gold particle: because of the high conductivity of diamond, A T g is slightly decreased by ~ 7 K, and the temperature 
in the diamond sphere A« 36 K is nearly uniform, approximately equal to the temperature at the centre of the 
same volume when replaced with water (see Figs. 16 z and [l6fc) . 

These numerical simulations can be compared with the analytical expressions derived above. According to Eq. 



FIG. 16. (a) Calculated temperature distribution in the xz plane with x g — 0 and a heating power Pheat = 36 fiW. (b) Same 
as (a) but with the diamond replaced with water, (c-e) Line cuts taken along the white dashed lines shown in (a,b). 
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FIG. 17. Temperature increase of the diamond nanosphere as a function of the position of the gold nanosphere, as obtained 
from the numerical simulation (red line) or from Eq. ( |25| . 


(21), the temperature of the gold particle in the presence of the diamond, A T g = 170 K, corresponds to an effective 
thermal conductivity of the system (glass + water + diamond nanosphere) of = 0.84 W-m 
to the conductivity of glass and water (0.8 W-m _1 -K _1 and 0.6 W-m _1 -K _1 , respectively. In Figs. 16 
the numerical simulations with the predictions of Eq. 
agreement. 


1 -K 1 2 3 4 5 6 7 8 . This is close 
e, we compare 


(22) using AT g = 170 K, which are found to be in excellent 
We can also run the simulation for various positions x g of the gold particle, taking into account the 


position-dependent heating power as given by Eq. (19). The resulting curve Td(x g ) is in very good agreement with 

In main-text Fig. 4e, the data were therefore fitted to Eq. (25) 
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the analytic formula (25), as illustrated in Fig. 
to extract the temperature of the gold particle A T g as well as the mean distance h g d . 
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